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We consider transport exponents associated with the dynamics of a wavepacket in a discrete 
one-dimensional quantum system and develop a general method for proving upper bounds for these 
exponents in terms of the norms of transfer matrices at complex energies. Using this method, we 
prove such upper bounds for the Fibonacci Hamiltonian. Together with the known lower bounds, 
this shows that these exponents are strictly between zero and one for sufficiently large coupling and 
the large coupling behavior follows a law predicted by Abe and Hiramoto. 
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I. INTRODUCTION 

The spreading of a quantum mechanical wavepacket 
has been the focus of intense research both in the physics 
and mathematics communities. The free case and, more 
generally, the case where the spectral measure of the ini- 
tial state has an absolutely continuous component are 
well understood. More recently, the case of point spec- 
trum has received a lot of attention and the question 
whether spectral localization leads to dynamical local- 
ization has been answered quite satisfactorily As a 
consequence, dynamical localization was shown in most 
systems for which spectral localization had been estab- 
lished earlier. The most prominent exception to this rule 
is given by the random dimer model |2(, which displays 
spectral localization and dynamical delocalization and, 
in fact, super-diffusive transport. 

Between these two extreme cases there are a lot of 
models for which intermediate spectral and dynamical 
behavior has been observed; for example, Bloch electrons 
in a magnetic field Q and one-dimensional quasicrystals 
0, H . The spreading of a wavepacket in such systems 
is much less understood. There are a number of works 
based on heuristics and numerics for the Harper and Fi- 
bonacci models devoted to this issue. As a consequence it 
is expected that transport in the Harper model is almost 
diffusive and transport in the Fibonacci model is anoma- 
lous in the sense that the transport exponents take values 
other than zero (which would be the case in a dynamically 
localized system), one-half (diffusive transport), and one 
(ballistic transport). In fact, these exponents decrease as 
the coupling constant is increased and their behavior for 
large values of the coupling constant A is expected to be 
of order 1/ log A 0,0. 

For the Harper model, there are no explicit rigorous 
bounds for the transport exponents. Bellissard et al. 
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show a lower bound in terms of the multifractal dimen- 
sions of the density of states, @, but they do not have 
lower bounds for these dimensions. On the other hand, 
there are some explicit rigorous results for the Fibonacci 
model. However, they are limited to bounding the trans- 
port exponents from below. The question of whether the 
expected upper bounds, and hence anomalous transport, 
hold has been open from a rigorous standpoint. In fact, 
there was no general method for proving upper bounds 
on transport exponents, and finding such a method (that 
is applicable to models of interest) is one of the most 
important problems in this field. 

In any event, quantum dynamics for models with sin- 
gular continuous spectra supported on Cantor sets with 
critical eigenfunctions is a rich subject and has been stud- 
ied by many authors, with most papers focusing on the 
Harper model and the Fibonacci model or generaliza- 
tions thereof. A nice discussion of "what determines the 
spreading of a wavepacket," particularly in these inter- 
mediate cases, may be found in [ji- 
lt is the purpose of this Letter to report on rigor- 
ous work concerning these issues. We have developed 
a general approach for proving upper bounds on trans- 
port exponents in terms of the norms of transfer matri- 
ces. These matrices are the standard tool, in one dimen- 
sion, to study spectral (and dynamical) properties of a 
given model. Essentially, it is always the ultimate goal to 
reduce a problem at hand (such as studying dynamical 
quantities) to properties of transfer matrices or, which is 
essentially equivalent, solutions of the associated differ- 
ence equation. The first rigorous results on lower bounds 
for transport exponents could only be established after 
Jitomirskaya and Last had found such a general corre- 
spondence 0] ■ See 0] for the first explicit lower bounds 
for transport exponents in the Fibonacci model and \H\ 
for the best bounds known to this date. We will describe 
our general method for proving upper bounds on trans- 
port exponents in Theorem ^ below. Secondly, we have 
applied this method to the Fibonacci model and proved 
upper bounds on transport exponents that are indeed 
of order 1/logA in the large coupling regime. See The- 
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orem [5] for the detailed statement. Together with the 
known lower bounds, this establishes the first rigorous 
anomalous transport result. 

The organizations is as follows: We recall the defini- 
tions of the transport exponents and some of the known 
rigorous lower bounds for them in Section [H] describe 
our general upper bound for these exponents in terms of 
transfer matrix norms in Section IlIII and discuss the ap- 
plication of this general result to the case of the Fibonacci 
model in Section HVI 



II. TRANSPORT EXPONENTS 

Consider a discrete one-dimensional Schrodinger oper- 
ator, 

[Hij>](n) = ip{n + 1) + 1/3(11 - 1) + V(n)<ip(n), (1) 

on l 2 ( Z) A number of recent papers (e.g., 0, 0, 0, 
0,0,0,0,01) were devoted to proving lower bounds 
on the spreading of an initially localized wavepacket, say 
tp = Si, under the dynamics governed by H, typically in 
situations where the spectral measure of Si with respect 
to H is purely singular and sometimes even pure point. 

A standard quantity that is considered to measure 
the spreading of the wavepacket is the following: For 
p > 0, define (|X|^)(T) = EnMM^T), where 
a(n,T) = f, e- 2t/T \(e- ltH Si,S n )\ 2 dt. Clearly, the 
faster )(T) grows, the faster e~ ltH Si spreads out, 

at least averaged in time. One typically wants to prove 
power-law bounds on (\X\ p Si )(T) and hence it is natural 
to define the following quantity: For p > 0, define the 
lower transport exponent /3j (p) by 

07 (p) = lid inf — — , — - 

and the upper transport exponent 0f (p) by 



obey < S + (a) < S~(a) < oo. These numbers control 
the power decaying tails of the wavepacket. In particular, 
the following critical exponents are of interest: 



0f{p) = limsup 



plogT 



Both functions 0f (p) are non-decreasing. Numerical 
studies suggest that for some models, they may exhibit 
non-trivial growth |3, [13, llSj . Such a multiscaling phe- 
nomenon was termed quantum intermittency by Guarneri 
and Mantica ^3- F° r the two main models of interest, 
numerics show that there is non-trivial growth for the 
Harper model, but no growth for the Fibonacci model 

Another way to describe the spreading of the wave- 
function, which turns out to capture the limiting behav- 
ior of 0f (p) for small and large values of p, respectively, 
is in terms of probabilities. We define time-averaged 
outside probabilities by P(N,T) — J2\ n \>N a ( n ^)- 
note S~(a) = — liminfT^oo los P [<^7^ 2 ' T ' > an d S + (a) — 
— limsupj^^ 1 ° sP [^ 7 7 2,T ^ for any a G [0,+oo]. They 



and 



= sup{a > : S ± (a) = 0} 



af = sup{a > : S ± (a) < oo}. 



We have that < aj < a~ < 1, < aj < a+ < 1, and 
also that aT < af, a~ < a+. 

One can interpret af as the (lower and upper) rates 
of propagation of the essential part of the wavepacket, 
and as the rates of propagation of the fastest (poly- 
nomially small) part of the wavepacket. In particular, 
if a > a+, then P(T a ,T) goes to faster than any in- 
verse power of T. Since a ballistic upper bound holds 
in our case (for any potential V), Theorem 4.1 in 0] 
yields lim p ^ P% (p) = af and lim^^ (if (p) = a± . In 
particular, since (3f (p) are nondecreasing, we have that 

(if (p) < ctf f° r every p > 0. 

When one wants to bound all these dynamical quan- 
tities for specific models, it is useful to connect them to 
the qualitative behavior of the solutions of the difference 
equation 

u(n + 1) + u(n - 1) + V(n)u(n) = zu(n) (2) 

since there are effective methods for studying these solu- 
tions. Presently, the known general results are limited to 
one-sided estimates of the transport exponents. Namely, 
a number of approaches to lower bounds on (if (p) have 
been found in recent years. The papers |lflll4lll5| derive 
such bounds in terms of the behavior of solutions to 
with real energies z, with a link furnished by Hausdorff- 
dimensional properties of spectral measures due to results 
of Guarneri, Combes, and Last |12|. In fact, what is 
proven in these papers (although not stated in this form) 
are lower bounds on . Therefore, the lower bounds for 
P~s t (p) obtained in this way are constant in p. There is 
also work by Guarneri and Schulz-Baldes 0] , who bound 
af from below in terms of the packing dimension of the 
spectral measure. 

Further developments of these ideas by Guarneri and 
Schulz-Baldes, Barbaroux et al., and Tcheremchantsev 
fig allowed these authors to obtain better lower bounds 
for pf (p) which are growing in p. These results elucidate 
the phenomenon of quantum intermittency. 

On the other hand, [ill develop a direct approach with- 
out an intermediate step. These papers use power-law 
upper bounds on solutions corresponding to energies from 
a set S to derive lower bounds for flj (p). The set S can 
even be very small. One already gets non-trivial bounds 
when S is not empty. If S is not negligible with respect 
to the spectral measure of S±, the bounds are stronger, 
but there are situations of interest (e.g., random poly- 
mer models |2(), where the spectral measure assigns zero 
weight to S. 
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III. TRANSFER MATRICES AND UPPER 
BOUNDS FOR TRANSPORT EXPONENTS 

It should be stressed that there were no general rigor- 
ous methods for bounding af~ , a J , or f3^ (p) non-trivially 
from above. In the present paper we propose a first gen- 
eral approach to proving upper bounds on af (which in 
turn bound af and P^ip) for all p > from above, as 
well) . This approach relates the dynamical quantities in- 
troduced above to the behavior of the solutions to the 
difference equation J5J) for complex energies z. To state 
this result, let us recall the reformulation of in terms 
of transfer matrices. These matrices are uniquely deter- 
mined by the requirement that 



< n + 1 )\ , J u(l) 

u{n) r^ B,z 'U(o) 



for every solution u of J2J). Consequently, 

{T(n, z) ■ ■ ■ T(l, z) n>l, 

Id n = 0, 

[rCn+l,*)]- 1 ---^,*)]- 1 n<-l, 

where 



T(m, z) 



z-V{m) -1 
1 



We have the following result: 

Theorem 1 Suppose H is given by {Tjl. where V is a 
bounded real-valued function, and K > 4 is such that 
o~(H) C \—K + 1, A" — 1]. Suppose that, for some C £ 
(0, oo) and a £ (0, 1), we have 

C H* E + *) II 2 ) 1 rfE = °( T " W ) 



IV. THE FIBONACCI HAMILTONIAN 

The Fibonacci Hamiltonian is an operator of the form 
JTJ, where the potential is given by V{n) — Xx\l-a,i) ( na 
mod 1) with a = (y/E — l)/2. This potential belongs to 
the more general class of Sturmian potentials, given by 
V(n) = Ax[i- Q .i) (na + 6 mod 1) with general irrational 
a £ (0, 1) and arbitrary 9 £ [0, 1). These sequences pro- 
vide standard models for one-dimensional quasicrystals. 
(See for the discovery of quasicrystals.) Early studies 
of the spectral properties of the Fibonacci model were 
done by Kohmoto et al. and Ostlund et al. Q|. It was 
suggested that the spectrum is always purely singular 
continuous and of zero Lebesgue measure. This was rig- 
orously established by Siito for the Fibonacci ca se l2J 
and by Bellissard et al. j^f and Damanik et al. |l4j in 
the general Sturmian case. Abe and Hiramoto studied 
the transport exponents for the Fibonacci model numer- 
ically @; see also Geisel et al. Q. They found that they 
are decreasing in A and their work suggests that 



const 
log A 



(3) 



as A — * oo. 

The general approaches to lower bounds for the trans- 
port exponent described above have all been applied to 
the Fibonacci Hamiltonian (and some Sturmian mod- 
els). The best lower bound for af was obtained by 



Killip et al. in It reads a l > 



where 



log 



v/TF 



C(A)+k+1/2> 

0.0126 and C(A), chosen so that one 



can prove a result like ^Li WHn, E)\\ 2 < CL 2<:< -^ +1 for 
energies in the spectrum of H (our definition differs from 

that of HI), obeys ((A) 



3 log Vs 



log I 



■ 5 + 1' 



(logA + 0(l)). This 



and 



shows in particular that a ; admits a lower bound of the 
type ©■ 

The best lower bound for a~ was found in |ll| , where it 
was shown that a~ > grrc^; • In terms of the exponents 

0$ (p), the best known lower bounds are (see |Tlj ) 



A" 



max 

3<-n<CT» 



$(n,£+£ 



dE = 0{T- m ) 



for every m > 1. Then off < a. In particular, (3g(p) < 
a for every p > 0. 

Remarks, (a) If the conditions of the theorem are sat- 
isfied for some sequence of times, — > oo, we get an 
upper bound for a~. 

(b) The proof of Theorem ^ is based on the well-known 

formula a{n,T) = ^ J \ ({H - E - i)" 1 ^ 6 n ) | 2 dE 
and the fact that resolvent decay is closely related to 
lower bounds on transfer matrix growth. Details will be 
given in [20| . 



P+2k 



/? 5 ~>) > { (P+ 1)(C(A) + K+1/2) 



P<2C(A) + 1, 
p>2C(A) + l. 



We also want to mention work on upper bounds for the 
slow part of the wavepacket by Killip et al. ^^|. More 
precisely, they showed that there exists a 5 £ (0, 1) such 
that for A large enough, P(C 2 T^ X \T) < 1-6. Here, 



6 log 



logS(A) 



and 



£(A) 



A- 4+ y/(X -4) 2 - 12 



(4) 



See [T^, Theorem 1.6. (i)] . However, this result does not 
say anything about the fast part of the wavepacket, and 
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in particular, no statement for any of the transport ex- 
ponents can be derived. 

With the help of Theorem^we can prove upper bounds 
for a+ for the Fibonacci model at sufficiently large cou- 
pling. These upper bounds show that ® is indeed true. 

The precise result is as follows: 

Theorem 2 Consider the Fibonacci Hamiltonian and 
assume that A > 8. Let ct(\) — 2 j°g^(A) > w hcre £(A) 
is as in (|4"|l. Then, a+ < a(X), and hence f3 + (p) < ct(\) 
for every p > 0. 

Remarks, (a) One can observe that a(X) < p(X). 

(b) Note that £(A) = A + 0(1) as A -> oo. Moreover, 

= 2 '"log 3 — ~ 0-876 and a(X) is a decreasing func- 
tion of A for A > 8. Thus, we establish anomalous trans- 
port for the Fibonacci Hamiltonian with coupling A > 8 



and confirm the asymptotic dependence of the transport 
exponents on the coupling constant A that was pre- 
dicted by Abe and Hiramoto. We emphasize again that 
this is the first model for which anomalous transport can 
be shown rigorously. 

(c) They key idea is to study the well-known trace map as 
a complex dynamical system and prove upper bounds on 
the imaginary width of the (complex) canonical approx- 
imants of the spectrum. Denote x n (z) — Tr$(i 7 ' n ,z), 
where F n is the n-th Fibonacci number and z G C, 
and o~ n = {z € C : |x„(z)| < 2}. Then we prove that 
a n C {z £ C : \Imz\ < C£(A)-"/ 2 }. Outside of the sets 
cr n U o~ n +i, the traces grow super-exponentially and this 
yields the lower bounds on the norms of transfer matri- 
ces that we need. The claimed upper bound on a+ then 
follows from Theorem ^ A detailed proof of Theorem [3 
may be found in |2Cj . 
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